Non-diffusive spin dynamics in a two-dimensional electron gas 
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We describe measurements of spin dynamics in the two-dimensional electron gas in GaAs/GaAlAs 
quantum wells. Optical techniques, including transient spin-grating spectroscopy, are used to probe 
the relaxation rates of spin polarization waves in the wavevector range from zero to 6 x 10 4 cm -1 . 
We find that the spin polarization lifetime is maximal at nonzero wavevector, in contrast with 
expectation based on ordinary spin diffusion, but in quantitative agreement with recent theories 
that treat diffusion in the presence of spin-orbit coupling. 
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Electronic systems with strong spin-orbit interaction 
(SOI) exhibit exotic effects that arise from the cou- 
pling of spin polarization and charge current, such as 
spin Hall cmv&ats^^^^ and current-induced spin 
polarization 8 . These effects involve manipulation of the 
electron spin via electric, rather than magnetic, fields, 
creating the potential for applications in areas from spin- 
tronics to quantum computing 9 . However, SOI is a 
double-edged sword, as it also has the undesired effect 
of causing decay of spin polarization, reflecting the non- 
conservation of the total spin operator, S, i.e. [S, H] ^ 0, 
where TL is any Hamiltonian that contains spin-orbit cou- 
pling. 

In the consideration of spin-orbit coupled systems, a 
great deal of attention has been focused on quantum wells 
or heterostructures fabricated in III-V semiconductors, 
where breaking of inversion symmetry allows coupling 
that is linear in momentum. For electrons propagating in 
[001] planes, the most general form of the linear coupling 
includes both Rashba 10 and Dresselhaus^ contributions: 



Hso = a(k y a x - k x a y ) + (3(k x a x - k y a y ), (1) 

where k x , y is the electron wavevector along the [10], [01] 
directions in the plane, and a and (3 are the strengths 
of the Rashba and Dresselhaus couplings. The spin- 
orbit terms generate an effective in-plane magnetic field, 
b so , whose direction depends on k. Spin nonconserva- 
tion takes the form of precession, at a rate governed by 
b so (k), during an electron's free flight between collisions. 
In the Dyakanov-Perel (DP) regime 1 ^, where the preces- 
sion angle during the free- flight time r is small, a spatially 
uniform spin polarization will relax exponentially to zero 
at the rate 1/t s cx |6 so | 2 t. Although this process relaxes 
the initial polarization state, spin memory is not entirely 
lost, even after the electron undergoes many collisions. 
The relationship between real-space trajectory and spin 
precession leads to a correlation between the electron's 
position and its spin. Such correlations are predicted to 
enhance the lifetime of certain spatially inhomogeneous 



spin polarization states beyond what would be expected 
for conventional spin diffusion 13 . 

Burkov et ali& considered Tt so with only Rashba cou- 
pling (/3 = 0), and showed that a helical spin density 
wave with wavevector ^/Thma/2 decays more slowly than 
a uniform (or q — 0) spin polarization, by a factor 32/7. 
This contrasts with ordinary spin diffusion, where the 
decay rate of a polarization wave increases monotoni- 
cally (cx q 2 ) with increasing wavevector. More recently, 
Bernevig et ulJ-L showed that the lifetime of the spin he- 
lix is enhanced further in the presence of both Rashba 
and Dresselhaus coupling, and diverges as a — > f3. The 
infinite lifetime, or "persistent spin helix" (PSH) state, 
is a manifestation of an exact SU (2) symmetry of Ti. so at 
the point a — (3. As a consequence of the ST/ (2) symme- 
try, a spiral spin polarization in the z, x+ plane (where z 
and x+ are the normal and [11] directions, respectively) 
with wavevector 4mev£+ is a conserved quantity. These 
predictions suggest the possibility of manipulating spin 
polarization through SOI, without necessarily compro- 
mising spin memory, by controlling a/ (3 with externally 
applied electric fields. 

Transient spin grating (TSG) experimentsi^i^ are par- 
ticularly well-suited to testing the theoretical predictions 
of Refsi 13 i 14 , as they directly probe the decay rate of 
nonuniform spin distributions. In the TSG technique, a 
spatially periodic polarization of the out-of-plane spin, 
S z , is created by a pair of interfering ultrashort laser 
pulses. By varying the relative angle of the two pump 
beams, we are able to vary the magnitude of the grating 
wavevector q in the range 0.44 — 5.3 x 10 4 cm -1 , cor- 
responding to wavelengths in the range from « 1 - 15 
microns. The subsequent time evolution of the spin po- 
larization wave amplitude, S^(t), is measured by coher- 
ent detection of a time-delayed probe beam that diffracts 
from the photoinjected "spin grating." 

In this work we report detailed TSG measurements 
on two samples, each consisting of 30 quantum wells of 
thickness 120 A. One was with <5-doped with Si in the 
barriers to create a high mobility (fi ~ 150, 000 cm 2 /V-s 
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at 5 K) two-dimensional electron gas (2DEG) with carrier 
density 7.8 xlO 11 cm~ 2 . The 5 layers are symmetrically 
located on either side of each well. The second sample is 
identical except that the mobility of the electrons was re- 
duced to ~ 3, 500 cm 2 /V-s by placing 83% of the dopant 
atoms in the wells, rather than in the barriers. 

The TSG measurements show that the lifetime of S* 
is maximal at nonzero \q\ in both samples, in contrast to 
expectations for simple diffusion but in agreement with 
the prediction of Burkov et al. 13 . The magnitude of the 
lifetime enhancement and its dependence on the direction 
of q point to the presence of both Rashba and Dressel- 
haus interactions in our samples. Quantitative agreement 
between the theory of Bcrnevig et ali^ and our exper- 
imental results suggests the feasibility of realizing and 
detecting a PSH in samples engineered to achieve a ~ (3. 
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FIG. 1: [Color Online] Decay of the out-of-plane component 
of spin polarization wave with wavevector in the [11] direction, 
for (a) the high-mobility sample with |g| = 0.58 x 10 4 cm -1 
and (b) the low-mobility sample with |<f| = 0.69 x 10 4 cm -1 , 
for several temperatures. The initial value of spin polarization 
is estimated at a few percent for all measurements. Decay 
curves are offset vertically for clarity. Solid symbols in panels 
(c) and (d) indicate the temperature dependence of the two 
time constants obtained from a double exponential fit to the 
decay curves. Solid lines are the lifetime of a spatially uniform 
spin polarization. 

Panels (a) and (b) of Fig. 1 show the decay of S q (t) 
for the high-/j, and low-^t samples, respectively, measured 
at the wavevectors where the decay rate is smallest. In 
both samples the decay of the grating is clearly not a sin- 



gle exponential at low T and crosses over to nearly single 
exponential as room temperature is approached. All the 
decay curves can be fitted to a double exponential form, 
a\ exp(— t/rx) + ai exp(— t/r^), with equal weighting fac- 
tors (a\ — 0,2) over almost the entire range of T. The 
only exception is the T < 25 K regime of the high-^i sam- 
ple, where the momentum relaxation rate, 1/r, becomes 
comparable to the spin precession frequency, f2. In this 
(f2r > 1) regime the initial decay is damped oscillatory 
rather than exponential 16 . The solid symbols in panels 
(c) and (d) are the time constants n and T2 extracted 
from the double-exponential fit, plotted as a function of 
T. For both samples the ratio between the fast and slow 
rates is at least a factor of ten at low T and gradually di- 
minishes as T approaches room temperature. The solid 
lines are the lifetimes of the uniform spin polarization, 
r s (T), as obtained from decay of transient Faraday rota- 
tion induced by a single, circularly polarized pump beam. 
For all temperatures, the q = polarization decays as a 
single exponential, whose lifetime lies between the two 
lifetimes observed at nonzero q. 

The dispersion of the decay rates with q\\ [11] is shown 
in Figs. 2(a) and 2(b), for the high and low \i samples, 
respectively. We present results for T = 50 K, a temper- 
ature at which both samples are in the fir < 1 regime, 
and yet the ratio t^/tx remains large (in a certain range 
of q). For both samples, the larger of the two time con- 
stants peaks sharply at |gj « 0.6 x 10 4 cm -1 . The life- 
time of the more rapidly decaying component decreases 
monotonically with increasing |<f|, consistent with simple 
diffusion. 

The spin-dynamical effects illustrated in Figs. 1 and 2, 
biexponential decay and non-monotonic dispersion, are 
in quantitative agreement with the theory of coupled 
charge and spin dynamics in the presence of Ti. so , which 
we describe briefly below. Assuming a single, isotropic r, 
Burkov et alM- and Bernevig et ali* derived a set of four 
equations that describes the coupling of electron density 
n q (t) and the three components of S q (t) brought about 
by the SOI. Along the [11] and [11] directions, the four 
equations separate into two coupled pairs. For q | [11], 
spin precesses in the z — x+ plane, leading to coupling of 
S x+ and S z . 

Solving the pair of equations that couple S z and S x+ 
yields two eigenmodcs and corresponding eigenfrequen- 
cies iuj\fi{l) = i/ t i,2(q)- The q || [11] eigenfrequencies 
are, 



sin 2(f) 



lor(q) 




± \l (1- sin 20) 2 + ^£(1 + sin 20), (2) 
% 

where <?o = m*v so /h is the reciprocal of the distance 
required for an electron to precess by 2-7T in the spin- 
orbit field and 70 = v 2 so k 2 F T is the typical DP (i.e., 
g^O) decay rate. The parameters <fi = tan _1 (a//3) and 
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FIG. 2: [Color Online] Decay lifetimes obtained from best 
fit of a double exponential to experimentally measured decay 
of spin polarization, as a function of wavevector q parallel 
to the [11] crystal axis. Data are shown for the (a) low-/i 
and (b) high-/i samples, at 50 K. The dashed lines are fits 
to a theory (described in the text) in which the two lifetimes 
correspond to the positive and negative spin helices illustrated 
in Fig. 2(a) inset. Fig. 2(b) inset: Comparison of spin 
polarization decay at q = 0, 0.58 x 10 4 cm -1 , and 2.01 x 10 4 
cm" . 



v so = \J a 2 + (3 2 reflect the relative and combined cou- 
pling strengths of the Rashba and Dresselhaus interac- 
tions, respectively. Finally, the spin diffusion coefficient 



\rj2 is given by 7o/2g, : 



D s 

The solutions for <f||[ll] are especially interesting at 
the SU{2) point where a — (3, and sin2c/> = 1. In this 
case the eigenvectors are (l,±i)/V5, corresponding to 
forward and backward spin spirals in the z,x + plane, 
for each value of q. A TSG experiment creates the ini- 
tial condition (1,0), which couples with equal strength to 
the two eigenvectors. The subsequent time evolution of 
the two eigenmodes has the form of a double exponential 
decay with equal weighting factors. At the "resonant" 
wavevector qo v8, the two decay rates are and I670. 
After decay of the unstable eigenmode, the PSH remains 
stable, despite the rapid electron scattering and spin pre- 
cession that are occurring. The eigenfrequenies for q\ \ [11] 
can be obtained from Eq. [21 as well, if we replace <fi by 
—<f). At the same SU(2) point where the PSH exists for 
q\\ [11], the spin dynamics for q\\ [11] obey simple diffusion. 



Several features of the physics at the SU{2) point re- 
main when a 7^ /3, but both are nonzero. The lifetime at 
the resonant wavevector is enhanced relative to the case 
when only one of two interactions is present, although it 
no longer diverges. The enhancement is expected to be 
stronger for <f|| [11] than q\\ [11]. The eigenvectors are still 
admixtures of S z and S x+ and the photoinjected wave of 
pure S z will again decay as a double exponential. How- 
ever, the weighting factors are g-dependent and the long- 
lived state will be an elliptical, rather than circular, spin 
helix. 

With this overview of the theory, we return to the ex- 
perimental data. The dotted lines in Fig. 2 show the 
best fits obtained by varying parameters in Eq. [21 The 
fits yield parameter values v so = 460 m/s, 4> = 0.2, and 
D s = 1000 cm 2 /s for the higher- [i sample and v so = 480 
m/s, 4> = 0.08, and D s = 105 cm 2 /s for the intentionally 
disordered sample. The above value of v so for the high-/! 
sample predicts a spin-precession frequency for ballistic 
electrons, f2 = v so kp = 0.10 THz, that is within 10% of 
the experimental valueit The theory applies equally well 
to both samples, despite a difference of about ten in their 
scattering rates. Note that smaller r leads to smaller spin 
relaxation rates, as the entire dispersion scales with the 
DP relaxation rate, 470. 

A surprising feature of the theorjii4 is the sensitivity 
of the dispersion curves to small admixtures of Rashba 
coupling into a pure Dresselhaus system, or vice versa. 
Our quantum wells are designed to be perfectly symmet- 
ric and therefore no Rashba term is expected. However, 
the maximum lifetime enhancement, ~ 7, for the higher 
/i sample is significantly larger than the factor 32/7 pre- 
dicted for a = 0. This discrepancy is accounted for by 
a relatively small admixture a ~ 0.2/3. It is possible 
that the presence of a small Rashba contribution can be 
traced differences in the interface on either side of the 
well, a known feature of the GaAs/GaAlAs systemi2*iS. 
It is interesting that the value a/ (3 ~ 0.08 is smaller in the 
lower-/i sample. The difference may reflect the tendency 
of the ionized impurities to attract electrons toward the 
barrier in case of <5-doping, as opposed to doping inside 
the well. 

A further prediction of the theory is that even small 
a/P generates a large anisotropy in the dispersion curves. 
To test this prediction we measured S z at 50 K in the 
higher-/i sample, for q oriented along the [11], [11], and 
[10] directions. The results shown in Fig. [3 demonstrate 
that the lifetime enhancement is strong for <t||[11] and 
weak for q\\ [11]. The line through the [11] curve is calcu- 
lated using Eq. [21 for the lifetimes and Eq. 22 of Ref. 14 
for the weighting factors of the double exponential decay. 
The parameters are the same as used to fit the disper- 
sion along [11]. The weighting factors and decay rates 
for q || [11] are calculated using the same parameters, 
provided that we replace <fi by — (p. The theory is clearly 
quite successful in providing a quantitative description of 
the spin relaxation dynamics seen in our experiment. 

The only feature of the data not captured by the the- 
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FIG. 3: [Color Online] Decay of spin polarization in the high- 
H sample as a function of time for three orientations of the 
grating wavevector (top-to-bottom: [11], [10], and [11]). The 
lines through the data are predictions of the theory described 
in the text, using the parameters that were obtained from 
fitting the dispersion curves. The anisotropy is smaller in the 
low- jj, sample, which is consistent with the conclusion based 
on the dispersion data. Inset: Fermi contour for a system of 
free electrons subject to the SOI described by Eq. Q showing 
the origin of the extremal orientations. Rashba coupling set 
to one-half of Dresselhaus (larger than in the systems under 
study here) for the purpose of illustration. 

ory described above is the gradual decrease in the Tijrx 
ratio with increasing T (shown in Fig. ^) that takes place 
for both samples. The characteristic scale of T clearly is 
not the spin-orbit splitting, which is ~ VfQo or about 1 
K for our samples. We speculate that the T-dependence 
of t-z/tx is a consequence of the cubic (in k) Dresselhaus 
coupling, HcD oc k x kya x — k y k 2 u y , which does not appear 
in Eq. ^ H c d is non-negligible when (k 2 ), the expecta- 
tion value of k 2 in the lowest subband, becomes compara- 
ble to kp. The cubic Dresselhaus coupling destroys the 
proportionality Q oc k, which ensures that the preces- 
sion angle in a free-flight between collisions depends only 
on the electron's displacement, and not on its velocity. 
In the presence of TicD the precession angle will depend 



on k, and the fractional spread in precession angle for 
a given Ak is ~ fci?Afc/(fc|). Assuming a thermally in- 
duced momentum distribution Ak ~ T/vp, the relative 
variation in precession angle is ~ T/Ei, the ratio of the 
temperature to the energy of the first excited subband. 
Our results are consistent with E\ (which is w 500 K 
for our samples) as the characteristic energy scale above 
which the correlations that generate the PSH are lost. 

Finally, we discuss directions for future research that 
are suggested by the conclusions of this work. The 
main implication of the theory and experiments described 
above is that SU(2) spin symmetry can be recovered, de- 
spite the presence of SOI, if the condition a = j3 can 
be satisfied. This can be accomplished by increasing the 
strength of the Rashba interaction, as compared with our 
samples, by designing asymmetric doping and/or barrier 
profiles, or by applying a gate electrode. Once achieved, 
the SU (2) symmetry is robust with respect to disorder 
and interactions^, and could be stabilized at room tem- 
perature by increasing E\. Schliemann et al. 20 previously 
discussed some attractive features of a system in which 
a — (3 for logic devices based on spin transport. They 
described a device structure in which spin is injected and 
detected at two points, and showed that if a = /3 the 
change in direction of the injected spin depends only on 
the vector displacement, a, of the two point contacts, 
even in the presence of strong scattering. Modulating 
the Rashba coupling via an external gate can then de- 
stroy the long-range spin correlation. The restriction to 
point contacts is actually unnecessary in the SU (2) sys- 
tem because the net spin precession depends only on the 
projection of a onto x+. For example, spin injected with 
arbitrary polarization from anywhere along a line contact 
defined by x + — will be detected with the same spin 
anywhere along a line contact at x+ = mr/2ma, where n 
is an integer. Thus, in principle, a planar spin logic de- 
vice operating at room temperature is feasible, if highly 
spin selective contacts to the 2DEG can be fabricated. 
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